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Abstract 

I consider the dilute monopole gas expansion of the three dimensional 
Yang-Mills-Higgs system in the symmetry broken phase. The functional 
determinants which occur in such an expansion are computed in the heat 
kernel approximation for an arbitrary SU{N) gauge group. Explicit ex- 
pressions for the gauge boson mass in the unbroken gauge sector and 
the string tension are obtained for the SU{2) gauge model and are eval- 
uated numerically. The results show a strong dependence on the ratio 
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1 Introduction 



The dynamics of three and four dimensional systems differ substantially. A char- 
acteristic difference is that the infrared modes play a more important role in the 
lower dimensional case. On the one hand this simplifies the perturbative investi- 
gations since the ultra-violet regime is less important and, by that, the separation 
of the non-physical, ultra-violet divergences becomes easier. On the other hand, 
however, the infrared modes may render the usual perturbation expansion diver- 
gent order by order in the infrared. In such a case we are forced to perform a 
partial resummation of the perturbation expansion. The resulting improved per- 
turbation expansion is based on certain quasi-particles with screened interactions 
and may be very different from the original perturbation expansion. 

These phenomena can nicely be studied by means of the three dimensional 
SU{N) Yang-Mills-Higgs system in the symmetry broken phase. The action of 
this model is given by 

S[Al $1 = J d'x (^^Ft^Ft^ + i(A$)"(A$)" + ^ - v^'Y^ , (1) 

where 

(A$)" = 9,$" + er''"4^$% (2) 

jafec g^j^g ^Yie structure constants of the SU{N) gauge group and the coupling 
constants and A have the dimension of an energy. Since the theory is super- 
renormalizable one only needs a single counterterm for the mass renormalization 
of the Higgs field But the perturbation expansion is order by order divergent 
in the infrared and thus renders the theory essentially nonperturbative. 

Fortunately, the action (|l]) possesses stable localized extrema, the static Yang- 
Mills-Higgs monopole configurations of the corresponding four dimensional gauge 
model 0, 1^, ^, 1^. The polarization of a dilute gas of such pseudo-particles 
screens the long range magnetic Coulomb forces so that the semiclassical expan- 
sion yields an infrared stable approximation . According to this expansion the 
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electric charges are confined by a linear potential and the gauge bosons acquire 
a non- vanishing mass. 

Originally, the reason for studying the SU (2) dilute monopole gas expansion 
was the fact that QED in the Georgi-Glashow model is a simple example of a 
compact U{1) gauge theory. Polyakov's work in 1976 had initiated an extensive 
study of the three dimensional compact U{1) gauge theory and in particular 
of a lattice version of this U{1) theory that is obtained by using the Villain 
approximation of the Wilson action . The primary goal of these lattice studies 
was and still is to reveal the phase structure of such theories by studying the 
behavior of quantities such as the string tension as a function of the inverse 
temperature P = 1/e^a where a is the lattice spacing. 

But there are various other important aspects of the three dimensional model 
(|l]) which require a detailed investigation of the dilute monopole gas expansion 
in the context of the model itself. One certainly is the fact that it resembles 
the effective model for the high temperature phase of QCD p|, ||, |T0|. Through 



dimensional reduction of the four dimensional gauge theory of strong interactions 
at finite temperature one arrives at a three dimensional gauge model in which 
the zero-component of the gauge field plays the role of the Higgs field. The 
Higgs sector is screened by the Debye mass whereas the gauge sector is infrared 
divergent. Therefore conventional perturbation expansion is inapplicable. To 
overcome this difficulty one has to perform a certain partial resummation of 
infinitely many Feynman diagrams. Such a resummation could be very similar 
to the dilute monopole gas expansion. 

Another important issue is the question of the true vacuum. In most com- 
putations of the effective potential only the trivial constant background field is 



considered. However it is well-known since the work of Coleman IT^, that 
nontrivial saddle point configurations of the underlying theory may alter these 
results fundamentally in the sense that what appears to be the vacuum state is 
actually unstable. For instance if we consider the three dimensional Yang-Mills- 
Higgs model in the symmetrical phase the one loop effective potential obtained 
with a constant background field shows no symmetry breaking effects. This may 
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change if we also include non-trivial saddle point configurations in the compu- 
tation. Such nontrivial saddle point configurations arise because of the explicit 
symmetry breaking source term that we commonly introduce in order to study 
spontaneous symmetry breaking. The form of these configurations is very sim- 
ilar to the 't Hooft-Polyakov monopoles and their contributions to the effective 
potential can be computed in the dilute gas expansion. However, in order to 
investigate such a possibility of spontaneous symmetry breaking we first must 
have a complete description of the dilute monopole gas. 

The present paper contains further details of the dilute monopole gas ex- 
pansion which was sketched in Ref. [Q. Its goal is to provide the formulae for 
the quantitative analysis of the expansion. This is achieved by an approximate 
computation of the 1-loop determinants which occur in the monopole partition 
function. We regularize the determinants with the help of the ^-function regular- 
ization and then evaluate the resulting expressions in the heat kernel approxima- 
tion [|1^, This enables us to present explicit expressions for the gauge boson 
mass in the unbroken U{1) gauge sector and the string tension in the SU{2) 
gauge model. 

The numerical evaluation of the formulas shows that the U{1) gauge boson 
mass and the string tension strongly depend on the ratio A/e^ = rn'^^^^^/m^Y 
where rriw is the heavy vector boson mass. In both limits, X/e"^ and A/e^ — >■ 
oo, the U{1) gauge boson mass and the string tension vanish but physics is 
different. We shall argue that in the Prasad-Sommerfield limit A/e^ the 
symmetry is restored whereas the limit A/e^ ^ oo corresponds to a non-confining 
Higgs phase with a massless U{1) gauge boson. 

Since there exist several misprints in the original literature and since the 
precise expressions which arise in the formulation of the dilute monopole gas 
expansion are important for our considerations we find it appropriate to present 
the dilute monopole gas expansion in greater detail. In Section 2 we rederive the 
one-monopole partition function for the SU{N) Yang-Mills-Higgs theory. The 
evaluation of the SU{N) one-loop determinants with the help of the heat kernel 
method is subject of Section 3. In Section 4 we give the precise expressions for 
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the SU (2) dilute monopole gas partition function, the gauge boson mass in the 
unbroken gauge sector and the string tension. Numerical results for the SU{2) 
gauge model are presented in Section 5 and a summary is given in Section 6. 



2 The SU{N) one- monopole partition function 

The general presupposition for the semi-classical expansion is the existence of 
non-trivial saddle point configurations of the underlying classical action. Since 
the action (|T]) of the three dimensional SU{N) Yang-Mills-Higgs system is just the 
energy of the static configurations of the corresponding four dimensional system, 
the static solitons of the four dimensional theory - the Yang-Mills-Higgs magnetic 
monopoles - are the saddle points or pseudo-particles of the three dimensional 
model. 

In order to obtain the one-monopole contribution to the functional integral, 

Z = I V[A1]V[^-] exp i-S[Al , (3) 

in the one-loop approximation one writes the fields as a sum of the 

classical one-monopole field A^ and a quantum fiuctuation a" (0°) and 

expands the action up to terms quadratic in the quantum fiuctuations. As usual 
the functional integration over the quantum fiuctuations requires some gauge- 
fixing. A suitable choice is the background-field gauge defined by the gauge-fixing 
action 

Sgi[A1, ^''] = \j d^x ipf (a\ - A^) + er^^%^ -^'yf ■ 
The corresponding Faddeev-Popov determinant is given by 

det [MM = det \Dt'Df + f-'^d jdc'h^c^c'^^ ^ 

which, in a one-loop computation can be approximated by 

det [A<Fp] ~ det [Sl^Df + f^<^d ^dc'b^c^c'A _ 
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If one performs the standard algebraic manipulations and employs the classical 
field equations 

(DiDi^Y = ^ (^"^ - v""^) , (7) 

as well as the Bianchi identity for the SU{N) structure constants one arrives at 
the following expression for the functional integral in the one-monopole sector 

= I V[at]V[r] det [Mfp] exp {-S^[Al - S^^^,[al r]) ■ (8) 
Here, Sin[A^, is the classical one-monopole action, 

S^[Al = J d?x {jPtiH + ^iMnOi^" + ^ (^"' - V^^f^ , (9) 

and 5'quadr[a", 0"] that part of the full action which is quadratic in the quantum 
fluctuations: 

S^.^A<.<t>''] = \jd'x{at, ~M^){f)^ (10) 

where 

j^ab ^ -Dl^Dl^ + ^ - v"'^^ (5"'' + A^''^* - . (n) 

Apart from the gauge zero modes which are eliminated by the gauge-fixing 
there are additional zero modes since the classical one-monopole solution violates 
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translational symmetry. The corresponding eigenf unctions of the quadratic form 
(|TT]) are given by 



Af = j d?x (f^^F^^ + {dST(dST) . (12) 

To ehminate these zero modes in the functional integral (D we again use the 
standard Faddeev-Popov method of insertion of unity. The appropriate decom- 
position of unity is 



1 = det 

kl 



X j dRp- / {alFt, + r{Di%r)^ (13) 

where R is the center of mass coordinate of the monopole. In the one-loop 
approximation the determinant in (|13D can easily be evaluated and gives 
Thus if we insert ([13|) in the functional integral @) we obtain 

Zi = j dn U^'^ det [Mfp] exp (-SJ^A", 

X j v[amr] n ^ / {atF[i + rimT)^ 

xexp(-5q„adr[<,0-]) . (14) 

Since no further zero modes exist we can now do the functional integral over 
the quantum fluctuations. If we use Pauli-Villars regularization the expression 
for the one-monopole contribution to the functional integral in the one-loop ap- 
proximation is given by 

Zi = JdR det [Mfp] (cfet [M])'^^^ exp [S^A^, l'^]) (15) 

where M is the Pauli-Villars regulator mass and the quadratic form defined 
in ( pUD and ([TT]). In case of the quantum fluctuation determinant we write det 



instead of det to indicate that the determinant has to be calculated with respect 
to the non-zero modes of Ai only. Furthermore it is implied that both functional 
determinants are regularized by the Pauli-Villars method and normalized by the 
corresponding 'free' determinants, det[A^5^p] and det[A^'^]. J^^pp and Ai^ are 
obtained from Aipp and Ai by setting A" = and $" = v"". 

Finally, it is useful to introduce the dimensionless fields f and — >■ 

vA°; as well as the coordinates Xi Xi/{ev), Ri Ri/{ev) and M Mev In 
terms of these dimensionless quantities the functional integral (|^) reads 

-1/2 



(16) 



where m\Y = ve is heavy vector boson mass. 



S[AM = '^ j d'x{-Ft^F^^ + -{mnm^ , (17) 



F^. = diA'^-djA1 + r'"A',A'^ , 



(A^)" = + /'^^'^^^^^ , (18) 



and 



^ab ^ -m^Df + ^ f*'^^ - l) d""^ + 4^"^* - . (20) 

2e ^ ■' 

With the help of the one-monopole result (|T^) one can now construct a formal 
expression for the partition function of a dilute gas of infinitely many monopoles 
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(see Section 4). However, since we are interested in an explicit expression for 
this partition function which then can be investigated numerically we are forced 
to work out the one-loop determinants. What we can say already at this point 
by only looking at the above formulas for the one-loop determinants is that the 
result will depend on the ratio A/e^, only. 



3 Calculation of the SU{N) determinants 

The method we are going to use to evaluate the determinants occurring in the 
expression for the one-monopole partition function is explained in detail in Refs. 
|l^ and IH^. It is based on the observation that the ratio of two determinants, 



det^ and deti3 with finite number of zero modes, n_4 and n^, can be written as 
/ det^\ dt 



t 



Tr (. 



(21) 



In general this expression has to be regularized. The most suited regularization 
scheme is the ^-function regularization defined by 

'det^\ _ 



In 



deti3/ 



reg 



M' 



2s 



Vis) Jo 



dtt 



s-l 



Tr (e 



-.-tA 



riA + riB 



(22) 



where Mis again the Pauli-Villars regulator mass. The functional trace can easily 
be evaluated with respect to the complete set of plane waves exp{—ipjXj). Then 
we only have to shift all differential operators dj occurring in A and B by ipj and 
integrate over d'^p. Consequently (^2D becomes 

det^' 



In 



deti3 



Ae. ^-0\ds [T{s) Jo 



d'^x 



(27r)' 



:tr e 



-tA(dj^dj+ipj) -tB{dj~^dj+ipj) 



(23) 
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where we used "tr" instead of "Tr" to indicate that the trace in the functional 
space has aheady been performed. 

Mostly it is rather hopeless to compute expression (^) exactly and one has 
to rely on a certain approximation. The so-called "heat kernel" approximation 
is based on a power series expansion in t of the integrand occurring in (PB)). 
For dimensional reason an equivalent method is first to expand the exponential 
functions in powers of the covariant derivative (or the classical background-field) 
and then to integrate over the momentum. The approximation consists of cutting 
off the power series expansion in t at a suitable order and replacing the upper 
limit of the t-integration by a finite limit to such that the integral is stationary 
at to- 



In Ref. [0 this method has been used to evaluate the determinants occurring 
in the functional integral of the pure gauge theory in the presence of an instanton 
background field. Thereby it turned out that already an approximation in which 
only the first two nontrivial terms of the expansion in t are taken into account 
agrees with the exact result within 3%. A possible explanation of the surprising 
accuracy of the approximation is that the first two nontrivial terms are already 
the third and fourth order terms in the t-expansion. 

Since the evaluation of the determinants becomes much more complicated 
for the Yang-Mills-Higgs system than for the pure gauge theory we shall restrict 
ourselves to the first two nontrivial terms in the t-expansion, too. However, it 
should be mentioned that in our case these contributions will result from the 
second and third order terms of the t-expansion. Therefore we cannot expect a 



similar accuracy like in the corresponding instanton computation of but we 
believe that the approximation should still be reasonable. 

Let us begin with the evaluation of the Faddeev-Popov determinant, i.e., 

(24) 

As described above we insert a complete set of plane waves and expand the 
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exponential exp(tA^Fp) in powers of the covariant derivatives. This leads us to 
Tre*^- = tr / d'x J ^e"^^* |l + + [$]^) t - Ap^p.D^D- 

-4piPj {{d^ + [$]2) D,D, + A {d^ + [$]') D, 
+D,D, (d^ + ^ + lQp,p^pkPiD,D,DkDi^^^ 

+0(^6)} 

= ^tr J d'x |r3/2 + [$]2rV2 + ^ (- AZ)2^, ^ A/^, Ai^, 

+5^[$]2 _ 2A[$]2a + + 3[$]^) t^/^ 

, (25) 

where, for simplicity, we already supressed all terms with an odd number of p^'s 
because they trivially vanish once the d^p-integration is performed. Furthermore 
we introduced the shorthand notation [X] which stands for 

^j^^ab ^ jacb^c ^26) 

and should not be confused with the commutator notation [Di,Dj\. 
The evaluation of (pSj) goes straightforward. One easily verifies that 

tr {d,D,D,D, - D^D'd) = ^tr ([A, D,] [A, A]) = , (27) 

and 

tr (A^ [$][$] + [$][$] - 2A [$][$] a) = 
tr ([ZXD$][$] + [$][DZ?$] + 2[D$][D$]) 

= -2N |^(A$)"(A$)" + ^ ($"^ - l) . (28) 

Substituting these expressions into ([25| ) and performing the remaining traces we 
arrive at 

(27r)3 J \ N 
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_^ ($-^y _ ^ ((i''^-|>''$'=)') tV2 + o(t=^/2)| . (29) 

The corresponding expression for the 'free' Faddeev-Popov matrix, A^pp is im- 
mediately obtained by setting A" = and = = t>"/t>. Thus 

Tr (e^FP* - e^FP*) = 



(2vr) 



(30) 

Since the operators Aipp and TWpp are positive definite and possess a con- 
tinuous spectrum we can employ the heat kernel approximation. Thereby we 
assume that the left-hand side of (|30|) is a rapidly decaying function of t which 
can be approximated by the first few terms of its expansion. Thus, to obtain 
an estimate for the t-integral in (p5[) we insert the expansion (|3^) into (^) and 
replace the infinite upper limit of the integral by a finite limit to- For to we choose 
that value at which the result as a function of to possesses an extremum. In this 
way we obtain 



/ det[-^pp] \ ^ liml"^ (^''^ f'dff^-'lcy f~''' 3 



+1/2 



(31) 



where 



"(2^ 



/v TT I r I \ X 

= W / "^'^ 14^5^5 + (a4)"(5.4)° + ^ ($°^ - 1) 
-|((*°')'-^)-3((''°"**)'-(''' 



(32) 
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and to, determined by the extremum condition a^pptg ^^"^ ~ ^Mff^o^'^ = is given 



by 



to = ^ . (33) 



IFP 



Note, that to must be positive. Otherwise we must take into account the next 
higher order in the t-expansion. In fact, we can easily convince ourselves that to 
indeed is always positive. Since < 1 the coefficient ajupp is clearly positive. 
As far as Pmpp is concerned the only negative term is the third one which, ac- 
cording to the field equation for $ will exactly be canceled by the second term. 
Consequently, both coefficients are always positive. 

Finally we have to perform the integral over dt and to take the limit s — ^ 
which leads to the following result: 

'det[->lFp] 



I ^„^r A^o^ 1 ~ 4:^(^MfpPmfp = 4:pMFpto^ ■ (34) 

reg 

Apart from the fact that the expressions become rather lengthy the compu- 
tation of the quantum fluctuation determinant, i.e.. 



reg 



(35) 

goes exactly the same way. Instead of ( [25| ) one finds 



Tre-^* = 

+ {(D' + [|)]2)5,, + 2[F,,] + V{^) + [<l]2) t 
-ApkPiDkDiiS,, + 1)^ + + myS,, + 2[F,,](D^ + 

+2Dk[F,,]Di + 2[F,,]DkDi + DkDi{D' - V($) + [^') 
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+ lQpkPlPmPnDkDiDmDni6ij + 1)- + 0(D^ 



TC 



3/2 



tr J d^x {4t-3/2 + (4[$]2 - V{$)) r 



1/2 



+ ^ (4 (D.DjDiD, - DiD^Di) - 12[Fij][Fi,] - 24[A$][A$] 
+4 - 2Di[^^Di + + 3[$]^) 

- (dV($) - 2Al^($)A + V{^)D'^ - 3F($)2) 
-3 ([$]V($) + [$]')) t^/' + 0{t^/^)} , 

(36) 

where, in order to simplify the expression at least a little bit, we used the short- 
hand notation V{^) for 



With the help of (|8D and 

tr + V{^)I? - 2AV($) a) 

(iV2 + 1) A 



e2 



(Ar2 + 1) A / p^^^ap.^^a ^ ^ /^a2 _ ^\ ^62 j (3^ 



g2 ^v--^ v--^ ■ 
we obtain from ( |5B[ ) 

Tre^* = 

_3/2 ^ 

- ('4jv$»' + A ($-2 _ n (iv^ - 1) + ^tA r"' 

+ ^ l^lOA^^^.^^. + 16A^(A$)''(A$)" - ^ ($'^2 _ ^62 

-(A^2 + 1)^(A$)"(A$)" + 24 {^''^y + 6N 



4e4 

+8 (^"^ - 1] ^^"^ + 12 ($"2 



)'^^ti/2 + o(t3/2)|. 



(39) 
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If we substitute this expression and the corresponding one for into ( p5D and 
remember that Ai has = 3 translational zero modes we arrive at 



(40) 



where 



IT 



4N + — (^N^ + 1U J (l - ^""^ 



(27r)3 



a2 1 \ /f.b2 



6(27r)3 7 \^ 



-(Ar2 + 1) — (Di$)"(A$)'' + 24 (^($"^) -1 

\2 

+ — (^''^ - l) ((iV^ - l) (^''^ - 3) + 4 (5$^^ + 3 
+6N (^[d''^''^^^'y - [d^^^^v'^y^^ . (41) 

The expression for ti follows from the extremum condition aMti ^^^+/5A4ii^^~3 = 
0. Since $^ < 1 the coefficient aM is always positive and therefore 



i}^' = " (42) 



3 



^, if/?^>0. 



1 /2 

Since t/ must be real and positive aj^ and must satisfy the inequality 
O-mI^m ^9/4. Otherwise the next higher order in the t-expansion must be taken 
into account. For N = 2, i.e., for the SU{2) gauge model, we have checked 
explicitly that the inequality holds. But so far we do not have a general proof 
that this inequality holds for monopole configurations of arbitrary SU{N) gauge 
groups. 

Under the assumption that the above inequality holds for arbitrary N we can 
integrate over dt and take the limit s 0. This then gives 

In iPjj^] - 37^ + 2aMh'^' - 2(3Mty' + 3 In {m\) . (43) 
Vdet[A^o]^,g ^ ^ 

where 75' ~ 0.5772... is the Euler constant. 
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If we insert (|3^ ) and (^) into ([T^) we obtain the following expression for the 
one-monopole partition function 



^1 



V e2 ; 



c / 



X exp 



1/2 



exp 



^e2/47r^ 



3/2 



/ c/R ;^A(A/e2) exp 



mw 



C{\/e^ 



where 



mw 



is the usual one-monopole mass integral and 



A{\/e^) = 2tx\ — 1 



3^ ■ 

- 



An J 



(44) 
(45) 

(46) 



Note, that in contrast to the corresponding result for the one-loop instanton 
partition function in four dimensions the Pauli-Villars regulator mass M 
drops out from the final expression. 

Having an explicit expression for one-monopole partition function at hand we 
can now study the dilute monopole gas. The simplest example of such a gas is 
the SU{2) monopole gas consisting of 't Hooft-Polyakov monopoles 



4 The SU{2) dilute monopole gas 

For the SU{2) gauge group and A/e^ > it has been shown that if the mean 
distance between two widely separated monopoles (|Ri — Rj| » 1) is large com- 
pared to (A/e^)^^ the superposition principle holds and the interaction between 



the monopoles is described by the Coulomb interaction |]16|. Thus in this case 
the action of a system of n monopoles with magnetic charges = ±1 can be 
written in the form 
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and the partition function Z for a gas of monopoles becomes 



E 



^1 



rfRfc exp 



E 



fc=i 



2e2 ^|R.-R„ 



(4J 



with (cf. (H)) 



1 



^1 



\ 3/2 

"""^ ^ A(A/e2)exp 



e2/47r 



C(A/e' 



(49) 



27r V(eV47r)^ 

Note, that in the Prasad-Sommerfield hmit A/e^ there appears a change 
in the interaction of monopoles |18], |T9[. The Higgs field becomes massless, 
and the attractive force associated with it becomes long range. This force cancels 
the repulsive magnetic force between like charges and doubles the attractive force 
between unlike charges. Consequently two Prasad-Sommerfield monopoles do not 
interact via the Coulomb force if they possess the same magnetic charge. 

The Coulomb gas of magnetic monopoles interacting with an external mag- 
netic charge p(x) is described by the partition function 



E 



^1 



J]^ c/Rfe exp 

fc=i 



E 



mi 



X 



e.p(>/AE.*-R.).(x)^ 



(50) 



where t] is the potential corresponding to the external magnetic charge, i.e., 
A?7(x) = 27rp(x). A path integral expression for Z[ri] is now easily obtained 
by rewriting the Coulomb interaction in (|50|) as an integral over an auxiliary field 
x(R). The resulting expression reads 

ld'x(^^iVix-v)y-^'cosx) (51) 



Z[r]] = J V[x]exp 
where u given by 



uj^ = Svr- 



-zi 



\ 5/2 

A(A/e ) exp 



mw 



eV4' 



TT 



C(A/e2 



(52) 



e- \_(eV4vr)^ 

is the gauge boson mass in the unbroken U(l) gauge sector in units of mw, that 
is, mu{i) = uJiTLw- 
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The Euclidean functional integral (^) describes a Debye plasma of monopoles 
and antimonopoles if the Debye radius given by 1/ 'mu{i) is large compared to the 
mean distance between the monopoles. If n is the density of monopoles the mean 
distance is proportional to n^^/^. The density of monopoles on the other hand is 
proportional to (em[/(i))^. Thus the condition for the Debye approximation is 

^ « 1 . (53) 

Since mu{i) is proportional to exp[— mi4/C(A/e^)/2(e^/47r)] it is clear that by 
taking (e^/47r) small the monopole gas can be made as dilute as one likes. The 
nonlinearities in functional integral (|5T|) become exponentially small, too, because 
^(7(1) /e^ is nothing but the dimensionless effective coupling of the monopole 
gas. However, because of the long-range magnetic Coulomb interaction between 
monopoles, the gas does not become noninteracting even for very small densities. 
This is a general feature of the three dimensional compact U{1) gauge model. 

The "area law" behavior follows from the partition function Z[ri\ if one con- 
siders an external monopole charge density pi(x) generated by an Wilson loop / 
in the xi, X2 plane. For a large Wilson loop / Pi{^) is approximately given by 
the monopole charge density which is generated by an classical external electric 
current loop ji = eijkdjHk with djHj = pi, i.e., 

p,(x) = 5'{xs)Qi{xi,X2) , (54) 

where Q(xi, X2) = 1 if xi and X2 are the coordinates of a point inside the loop / 
and zero otherwise. 

In the semi-classical approximation Z[rii\ is approximated by 

Z[„A « /d|x1 exp [-j^jA (i(V(x. - 

where Xci is a solution of 

V^Xci = V^r]; + u"^ sin Xci = Svrpi + uj"^ sin Xci • (56) 
An approximate solution of (|S§) with pi given by (Bl) is 



^ I /V7/., _ ^ COSXcl 



(55) 



49i(xi, X2) arctan (e '^^3), if xs > 0, 
Xci = { (57) 
49i(xi, X2) arctan (e'^^'^) , if X3 < 0. 
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If we insert this solution into ( |55| ) we finally arrive at 



(58) 



where Si is the enclosed by the loop / and 

2e2 



cr 



2e^ 



(59) 



the string tension. This is Polyakov's formula (corrected for misprints). 

To summarize, the polarization of a dilute monopole gas of 't Hooft-Polyakov 
monopoles screens the long range magnetic Coulomb force so that an infrared 
stable semiclassical expansion becomes possible. There are no massless particles 
in the theory and non-singlet states are confined by a linearly rising potential. 
The U{1) gauge boson mass and the string tension of the dilute gas are given by 



ijj 



mw A 



_eV47r)' e2 



and 



where 



UJ 



mw A 



eV47r) 



mw A 



eV47r) 



mw 
(e2/47r) 



5/2 



A{X/e^) exp 



mw 



C{X/e') 



(60) 



(61) 



(62) 



5 Numerical results for the SU(2) gauge model 

The monopole configuration of the SU (2) gauge model is of the form [|l| 

A{r)a'' = ax f ^ , 



$'^(f)a'^ = f-a^, (63) 
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where a"" are the Pauh matrices. The functions h{r) and g{r) are obtained by 
minimizing the action (p!7D, that is the monopole mass integral 

C{\/e^) = rdr\ {rh\r) + h{r)f + 2h\r) { 1 + 



+ '^g'\r) + g\r) (1 + rh{r)f + (/(r) - l)'^ . 

(64) 

The corresponding differential equations are 

r^/i"(r) + 2rh\r) - 2h{r) - rg^{r) (1 + rh{r)) - 3rh^{r) - r^h^{r) = , 



r'^g"{r) + 2rg'{r) — 25f(r) + -;^r'^g{r) — Arh{r)g{r) 

A 



A_ 
2e"2' 



-2r^h^{r)g{r) ^r2/(r) = 



(65) 



which apart from the Prasad- Sommerfield limit (A/e^ 0) can only be 
solved numerically. 

To determine the functions h{r) and g{r) we use instead of (^) the cor- 
responding system of coupled nonlinear integral equations derived by Bais and 
Primack which can be solved numerically by simple iteration. In terms of 
h{r) and /(r) = 1 — g{r) these integral equations are given by 



h{r) = ho{r)+ Gh{r,r'){Ah{r') + B[ho{r')])dr\ 
Jo 

POO 

fir) = goir)+ Gf{r,r'){Af{r') + B[foir')])dr' , (66) 







where 

/ cosh(r) _ sinh(r) \ e''^' jl+l/r') r' > r 

( coshjr') _ sinh(r') \ e~-{l+l/r) jf / < ^ 



Ghir,r') 
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G/(r,r') 



I V 



V 

^ coshC^AT^r') _ smh{./x/^r') \ e'V^' (l+l/ ( y/y^^r) ) 



if r' < r, 
(67) 

are the Green's functions of the homogeneous equations corresponding to (0), 
Ah{r) = -r{l+rh{r)){f^{r) -2f{r)) -3rh\r) -r^h^r) -r, 

Af{r) = 2 + 2r/i(r)(2 + r/i(r))(l-/(r)) + A^V'(r)(3-/(r)), (68) 
are the nonhnear sources and 

B[ho{r)] = r^h'^{r) + 2rh'Q{r)-{2 + r^)ho{r), 



B[fo{r)] = rV^(r) + 2r/^(r)- (^2 + -r^J/o(r). 
The functions hQ{r) and /o('") must satisfy the boundary conditions 

f(r) — > 1 - Cf(X/ey + 0(r^), hir) — > Ch{Xle^)r + ©(r^), 



(69) 



but can otherwise arbitrarily be chosen. As in [ET] we have used 



and /o(r) = 1 



(71) 



r2 + 6(A/e2) ^ ^ (r2 + a(A/e2))i/2 ^ 

where a(A/e^) and 6(A/e^) are to be determined by minimizing the monopole 
mass integral C{\/e^). The functionals -B[/io(^)] and i?[/o(r)] then take the form 

26(A/e2)r(3r2 -6(A/e2)) 



S[/io(r)] 



(r2 + 6(A/e2))3 
a(A/e2)r(r2 - 2a) A ^2 



(2 + r2)/io(r) 



(r2 + a(A/e2))5/2 



- -r2 + 2 /o(r). 



(72) 
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For a detailed derivation of the coupled set of integral equations ( |66D we refer the 
reader to ref. [ pT[] . 

To study the A/e^-dependence of the U{1) gauge boson mass and the string 
tension a we have calculated iterative solutions of the system ( pB| ) for various 
values of A/e^. The results for the mass integral C(A/e^) and A(A/e^) which 
together determine the total A/e^-dependence of m;7(i) and a are plotted in Fig. 
1 and Fig. 2. In addition we have listed a few illustrative values in Table 1. 



A/e2 


0.1 


0.5 


2.0 


10.0 


100.0 




1.106 


1.119 


1.291 


1.433 


1.617 


A{X/e') 


3.7-10-19 


9.2 - 10-3 


6.576 


1.74-10-1 


1.14-10-1° 



Table I. C{X/e^) and ^(A/e^) for different values of A/e^ 



In contrast to the mass integral C(A/e^) which is a rather slowly increasing 
function of A/e^ between 1.0 (A/e^ = 0) and 1.787 (A/e^ oo) the function 
A(A/e^) varies substantially. This strong A/e^ dependence can formally be un- 
derstood if one compares the various terms which contribute to mass integral 
C(A/e^) and the coefficients aMpp, Pmpp, oim and [3m in the heat kernel expan- 
sion. Expressions that do appear in these coefficients but not in the mass integral 
are for example / ^^^(l — $2) and / d^x{{(^'^Y ~ 1)- contrast to the similar 
term (A/e^) / d?x{(^'^ — 1)^ which appears in the mass integral the size of these 
expressions depends very much on the shape of Higgs field which is determined 
by the ratio A/e^. 

At A/e^ ~ 2.2 the function A{\/e^) has a maximum. Since on the other hand 
A/e^ = 'mHiggs/''^H/ this result is equivalent with the statement that A{\/e^) is 
largest if the square of the Higgs mass is roughly twice as large as the square of 
the heavy vector boson mass. For A/e^ > 2.2 A{\/e^) decreases again and seems 
to vanish in the limit A/e^ — > oo. 

The behavior of A{\/e^) for small and large values of A/e^ becomes more 
transparent in Fig. 3 where the logarithm of A{\/e^) is shown. For curiosity we 
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have also tried to fit the curve for small and large values of A/e^ and arrived at 
the following result. For small A/e^ the function A{\/e^) seems to starts out as 

A(A/e2)^ 125.66(A/e2)i/=^exp(^-^j , (73) 

and the behavior for larger A/e^ may be described by 

expf-1.21(A/e2)2/3) 
A{\/e') ^ 100.56 ^ ^^^^;^J3 ^ . (74) 

However, at the present stage we should not pay too much attention to the precise 
values of the constants in (|73D and (|7^ since we do not know yet how accurate 
our approximate calculation of the determinants really is. 



Since and a are proportional to \J A{\/e'^) the above results show that 

both quantities vanish in the limits A/e^ ^ and A/e^ oo. However, physics 
is very different in these limits. Let us consider the case A/e^ oo first. For 
large values of A/e^ the attractive force associated with the Higgs field is clearly 
short range since the potential term in the classical action (0) forces the Higgs 
field to unity, i.e., the Higgs field rises from at r = to 1 within a distance of 
the order of (A/e^)~^. Therefore the interaction between the monopoles should 
indeed be very well described by the magnetic Coulomb force. The result that 
m[/(i) and a become small as we take A/e^ large becomes plausible if we consider 
the energy of an assembly of monopoles [|TB|. Classically, this energy is given by 
the minimum of the three dimensional action ([l|) 

+^|$|2(A$)"(A$)" + ^ (|$r - v^f^ . (75) 

subject to the constraint that the normalized Higgs field = |$|^^$" satisfies 
the required asymptotic properties of magnetic flux and nontrivial homotopy, i.e., 

lei,kr''d,^''dj^'dk^^ = 4rrJ2 QnHR - R„) . (76) 

2 n 

Since according to ( [76D di(f) possesses a 1/r singularity at the position of each 
monopole the quantity mina;[(Dj$)"(Dj$)"] should increase as the monopole den- 
sity increases. However, it is easy to see that if we take A large the potential term 
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forces ~ and it will be energetically favorable for mina;[(Dj$)"(Dj$)'^] to 
be small. Thus for large values of A the density of monopoles is small. 

If we go to very small values of A/e^ the situation is more complicated. First 
we must check that the superposition principle is still accurate, i.e., the attrac- 
tive force associated with the Higgs field can be neglected. This is the case if 
the distance between the monopoles in the gas is large compared to (A/e^)~^. 
Since the mean distance between the monopoles is proportional to (em[7(i))~^/'^ 
(cf. previous section) and according to the above results is proportional to 

(A/e^)^/'^ exp[— const. /(A/e^)] the superposition principle indeed seems justified 
even for very small values of A/e^. To understand the drastic decrease of m[/(i) 
and cr for small values of A/e^ let us again consider the classical energy of an as- 
sembly of monopoles (|75D subject to the constraint (^). We assume the presence 
of an extremely small monopole density and start to increase it. As the monopole 
density increases the quantity mina;[(Dj$)"(Dj$)"] increases. But when this term 
becomes larger than Af^ it will be energetically favorable for |$| to vanish which 
now is possible since the potential term in ([73| ) does not constrain |$| for small 
values of A. But once |$| vanishes it will be energetically favorable for the gauge 
fields to vanish, too. 

Indeed a detailed study of the A/e^ dependence of the expressions aj^pp , Pmfp ; 
and (3^ shows that in the limit A/e^ —>■ oo those terms with the highest power 
in the Higgs self-coupling A/e^ become dominant. These are just the terms which 
result from the Higgs potential. Since they appear with a negative sign in the 
exponential function (|46|) A(A/e^) becomes small as A/e^ becomes very large. In 
the limit A/e^ — > 0, on the other hand, it is not so that y4(A/e^) vanishes because 
the Higgs coupling becomes small but because the core of size l/(A/e^) outside 
of which the Higgs field approaches its asymptotic form becomes large. 

Thus based on the assumption that for small values of the coupling e^/47r the 
dilute monopole gas partition function is a reasonable approximation to the path 
integral of the Yang-Mills-Higgs theory the quantum theory undergoes a smooth 
transition to a non-confining Higgs phase with a massless U{1) gauge boson in 
the limit A/e^ — > oo and a rapid but still smooth transition to a non-confining 
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symmetrical phase in the hmit A/e^ — > 0. However, for any finite value of A/e^ 
the quantum theory is confining and all fields massive. 

To illustrate both dependencies namely, that on e^/47r and A/e^, the ^7(1) 
gauge boson mass is shown in Fig. 4a and Fig. 4b as a function of e^/47r for 
various values of A/e^. (The corresponding graphs for the string tension look 
very similar.) Since m[/(i) is proportional to exp[— mwC(A/e^)/2(e^/47r)] we see 
a rapid decrease for very small values of e^/47r. But we also see that the curves 
intersect. The reason for this is that the A/e^-dependence of m{7(i) enters in two 
different ways, namely on the one hand via A{\/e'^) and on the other hand via 
exp[— miyC'(A/e^)/(e^/47r)]. Although the A/e^ dependence of the mass integral 
is rather weak as compared to A{\/e^) it can play a certain role for sufficiently 
small values of e^/47r. To see this let us compare mu{i) for two different values 
of A/e^, say ci and C2 with ci < C2 < 2.2. For A/e^ < 2.2 A{\/e^) is an increas- 
ing function of A/e^ whereas exp(— mvi/C(A/e^)/(e^/47r) is always a decreasing 
function. Thus if e^/An becomes sufficiently small mu{i){c\) > m[/(i)(c2) but 
otherwise m[/(i)(ci) < mu{i){c2)- 



6 Summary 

The numerical results presented in the previous section show that it is indeed 
worth studying the dilute monopole gas expansion of the Yang-Mills-Higgs model 
in more detail, especially, if one views the expansion as an approximation of the 
Euclidean path integral of this model. 

The dilute monopole gas expansion of the Yang-Mills-Higgs model depends 
on two parameters, e^/Air and A/e^. The e^/47r-dependence becomes already 
transparent in the formal expression for the monopole partition function whereas 
the A/e^-dependence can only be seen when both the classical one-monopole 
mass integral and the corresponding functional determinants have been computed 
explicitly. 

In this paper the determinants have been evaluated with the help of the heat 
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kernel expansion, an approximation scheme which has been proven successful in 
a corresponding calculation for instantons. The result of this calculation is an 
analytical expression which we denoted by A{\/e^) and which is valid for arbi- 
trary SU (N) gauge groups. However, since apart from the Prasad-Sommerfield 
limit (A/e^ 0) the monopole solutions are not known analytically the further 
evaluation of A{X/e^) has to be done numerically. 

As an example we have considered the SU (2) gauge group. Although the mass 
integral C(A/e^) is only a rather weakly increasing function between 1.0 and 1.787 
its A/e^ dependence plays a certain role for small values of the couphng e'^/An. 
The much stronger A/e^ dependent quantity is A{X/e'^) which decreases for small 
and large values of A/e^ and has a maximum around A/e^ ~ 2.2. For any finite 
value of A/e^ the dilute monopole gas is confining. However, the string tension 
and the generated U{1) gauge boson mass are in general very small. In both 
limits, A/e^ oo and A/e^ — * 0, the system becomes non-confining but for 
physically different reasons. For A/e^ — >■ oo we see a smooth transition from the 
confining phase to the non-confining Higgs phase and for A/e^ ^ a rapid but 
still smooth transition from the confining phase to the non-confining symmetrical 
phase. 

Although we have seen that the superposition principle seems accurate even 
for very small values of A/e^ one would prefer a direct check. However, this 
would require to take into account the attractive force associated with the Higgs 
field which is by far nontrivial. What can certainly be done is to examine the 

accuracy of the heat kernel approximation used in this paper. To this extent one 
must compute the next order, i.e., the fourth order in the t-cxpansion and see 
whether these higher corrections alter the present results significantly. 
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Figure Captions 

Figure 1: The monopole mass integral as a function of A/e^ 

Figure 2: Tlie A/e^-dependence wtiicli results from the computation of the one- 
loop corrections to the monopole partition function. 

Figure 3: The logarithm of y4(A/e^). A fit for small and large values of A/e^ 
leads to an asymptotic behavior as given in (^) and (|7^ . 

Figure 4a: The logarithm of the U{1) gauge boson mass in units of Tfi^ as a 
function of e'^/An in units of mw for various values of A/e^; 

: A/e^ = 0.1, : A/e^ = 0.4, : A/e^ = 2, 

: A/e^ = 15, : A/e^ = 100. 

Figure 4b: See Fig. 4a. 
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